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The expression of the free energy density of a classical crystalhne system as a gradient expansion 
{2J[J' in terms of a set of order parameters is developed using classical density functional theory. The goal 

I here is to extend and complete an earlier derivation by Lowen et al (Europhys. Lett. 9, 791, 1989). 

, The limitations of the resulting expressions are also discussed including the boundary conditions 

' needed for finite systems and the fact that the results cannot, at present, be used to take into 

' account elastic relaxation. 
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I. INTRODUCTION 



Classical Density Functional Theory (DFT) can be used to determine the thermodynamic and structural properties 
' of classical systems based on only the interaction potential of the atoms and the external fields (including those 
of the walls if any) they are subject to[l]. It is based on the fact that the Helmholtz free energy is a unique 
^ ; functional of the local density and the goal of practical formulations of DFT is to provide an approximate form for 
: that functional relation that can be used to calculate the properties of real systems under a variety of thermodynamic 
^ ' and geometric conditions [1], [2]. For inhomogeneous systems, containing e.g. surfaces or bulk interfaces, even the 
I simplest approximate functionals can be analytically complex to work with so that from the earliest days of the 
■ O development of DFT in the 1970's it has been common in such circumstances to use gradient expansions in the order 
^ \ parameters in order to obtain simplified forms of the theory[l]. In fact, two types of quite different gradient expansions 
are currently in use. The first, described by Evans[l] and refined by Haymet and Oxtoby[3],[4], yields equations for 
the local density which vary with the microscopic structure of the system: for example, when applied to a solid, 
that some terms in the equations will exhibit variations on the length scale of the lattice parameter. The second 
type of gradient expansion is that of Lowen, Beier and Wagner [5, 6] which yields a typical Ginzburg-Landau model 
free energy. While based on similar ideas to the first method, this results in a true long-wavelength approximation 
' in which microscopic details, such as the underlying lattice when applied to a solid, are subsumed into the various 
\ coefficients of the gradient expansion. The resulting simplified free energy functional therefore takes the form of a 
OO ■ typical phenomenological Landau theory, but with coefficients that can be calculated from microscopic information. 
■ This theory has been used to discuss surface melting, the nucleation of freezing in a bulk fluid [7], [8] and solid-solid 
. phase transitions [9]. 

The purpose of this paper is twofold. The first goal is to review the two types of gradient expansions in order to 
clarify the differences between them. The need to do so is motivated by certain criticisms that have been made in 
the literature[10] about the inadequacy of the gradient expansions: as shown below, these seem likely to be relevant 
C ' to the older class but not to the Ginzburg-Landau model. The second goal is to clarify, extend and complete the 
. derivation of the Ginzburg-Landau models. The derivation currently in the literature [5], [6] was developed only for 
^ ' planar interfaces and is also incomplete. Furthermore, issues such as the appropriate boundary conditions for finite 
O ' systems have not been addressed. 

O In the next Section, the required elements of classical DFT are reviewed and the older form of the gradient expansion 

J> ' is discussed. Section III presents in detail a derivation of the Ginzburg-Landau model of Lowen et al. Finally, in the 
last Section the differences between the two types of expansion are discussed, as are the advantages and limitations 
of the models. 



II. DENSITY FUNCTIONAL THEORY 



Consider a system of N atoms with positions and momenta qi and pi respectively and mass m^. The atoms 
interact via a pair potential v (T*) and are also subject to the effect of an external field (p (^) which includes any 
one-body forces such as gravity or applied electrical fields as well as the effects of any walls. The local number density 
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IS 

N 



r(r)=^(5(r^-t.) (1) 



i=l 

and its average in a grand-canonical ensemble is denoted 

P(^) = (P(^))- (2) 

The grand distribution is 

where h is Planck's constant, /? = l/fcsT is the inverse temperature, is the chemical potential and the Hamiltonian 
is 

^ = fl^P' + 12^(^ii) + j dl'p{-f)cl>{-f). (4) 

The grand partition function is 

- = E / ^t^d exp {-PH + MiV) (5) 

and it is related to the grand potential, f2, as 

n = -kBTlnE. (6) 

The key results of DFT[1],[2] are that (a) there is a one-to-one correspondence between fields (j){~r^) and average 
density distributions p (V) leading to the fact that (b) there exists a unique functional of the density F [p] with the 
property that the functional 

([n] , [(!>]) = F[n]- n j dPr n{'r) + j dPr n {'r) cf) ( r") (7) 
is minimized, for a fixed field, by the equilibrium density 

(n = p) = (8) 



6n 



where the notation indicates that the variation is performed at fixed external field (j) as well as fixed chemical potential 
and temperature [1] , [2] . The functional F [n] is commonly written as F [n] = Fid [n]+Fex [n] with the ideal contribution 
given by 

pFia [n] = Jdr {n (7*) In (A^n (7^)) - n (7*)) , (9) 

and where the excess contribution is in general unknown. However, it is related to the m-body direct correlation 
functions by 

'-^ = - SniT',)...Snir^y ^''^ 
Combining eqs.(7)-(10), the Euler-Lagrangc equation can be written as 

In (A^n ("r )) - a ( r*) - H- /3(/) ( r*) = 0. (11) 
Finally, in terms of these quantities, the Helmholtz free energy A is 

A = F\p] + j dPr p{-r)(j){-r). (12) 
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All practical uses of DFT arc based on approximate functionals F^x [n] motivated by various physical reasons. 
Given such an approximation, one then has two choices as to how to proceed. The first is to use the functional to 
derive the one-body direct correlation function from eq.(lO) and then to derive the density from eq.(ll). (Note that 
in bulk systems, the only external force of interest is often that associated with the walls which contain the system. 
As this volume goes to infinity, it is expected that interior points are unafi^ected by the wall-force so that in these 
so-called "self-sustaining" systems, one can take (3(j>{~r) = 0.) The second approach is to parameterize the density 
and to then determine the parameters via eq(8). To be specific, let the parameterized density be n (V) = p(T^;r) 
where F = {Ta} represents a family of parameters. In a bulk system, the parameters are constants and eq.(8) implies 
that 

dF dN f ,^ dprr;r) 
If the parameters depended on position, the result would be 



SF 



Sn{r) 



^ M +4'[r) = 0. 14) 

p(T^;r) dTair) 5r„ ( r ) 9r„ ( r ) 



One trivial example of such a parameterization is a bulk liquid in which case one takes p{l^;T) = po for some 
constant pQ. Then, since the field is taken to be zero in the interior, the Euler-Lagrange equations simply give the 
usual thermodynamic relation 

dF 

^-,y = o. (15) 

A less trivial example is that of a bulk crystalline solid which is typically parameterized by a sum of Gaussians giving 



p (y ; r) = ^ (^) exp (-a - R?) 



(16) 



where {Ri (po)}^!^! ^^'^ the Bravais lattice vectors within the volume V and are, as indicated, functions of the average 
density po- Then, the parameters are F = {a,po} and since the field is taken to be zero in the interior, the Euler- 
Lagrange equations reduce to 



dF 
da 

dF 

dpo 



(17) 
-liV = 0. 



A solid-liquid interfacial system could be described by taking, e.g., a to be spatially dependent and by using eq.(14). 
Note that this requires that the free energy functional F [n] be known not only for bulk systems but for inhomogeneous 
systems as well. Given this functional, the choice of whether to use a parameterized density or to directly solve eq.(ll) 

is a matter of convenience. 



A. Gradient Expansion 



The classical gradient expansion of DFT was described by Evans[l] and given a more modern derivation by Haymet 
and Oxtoby[3],[4]. Here, the latter derivation is reviewed so as to contrast with the Ginzburg-Landau theory discussed 
below. 

Equation(lO) can be functionally integrated in density space. Define a parameterized nx (T*) = no + 
A (ni (T*) — no)where rii ("?*) = p (T*; F (7*)) is the state of interest and no corresponds to any convenient uniform 
(i.e. liquid) state. Then two integrations and adding and subtracting the same expression for m {'r) = ni gives 

f3Fex[ni]- /3Fea:{ni) = - d~rid~r2 / / C2 ("r^i, 7*2; nA') ("i (^i) - "o) (rii (^2) - no) dAdA'. 

J Jo Jo 

+ 61^1(0^2 / C2{l^i,'r2;nx'){ni-no){ni-no)dXdX'. (18) 

J Jo Jo 
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Adding in the ideal contribution gives 

pF [m] - l3F (ni) = y rfr*i (m (yOlnni (r*i) -nilnni) (19) 

- / rf7*irfr*2 / / C2 ( r*i,7*2;nv) (m ( - no) (ni (7*2) - no)dA(iA'. 
J Jo Jo 

+ I (rf\dPr2 I / C2 (7*1,7*2;^^') (m - no) (ni - no)dArfA', 

J Jo Jo 

and it should be noted that all of these relations are exact. 

Haymet and Oxtoby then introduce the following approximations. First, the two-body dcf is expanded about that 
of a uniform system as C2 (7*i, ~r2', ny) = C2 (ri2; no) + ...and higher order terms dropped giving 

pF [n] - pF (ni) y (i r^i (m ( T^i) Inm ( r^i) - ni Inni) (20) 

y" C2 (ri2; no) (ni (T^i) - ni) (m (7*2) + ni) d"r id"r*2 

Introducing the local free energy difference for a uniform system 

A/3/(7*i;ro) = pi(yi;ro)lnpi(yi;ro)-pilnpi (21) 

y C2 (ri2; no) (pi ( r*i; Tq) + Pi) (pi (7*2; Tq) - p^) d7*2 

the last result can be written more suggestively as 

/3F[n]-/JF(ni) ~ j i Ap f v,T i)) 

y C2 (ri2; no) (pi (7*1; To) + pj (pi (7*2; To) - pj dtid'r^ (22) 

so that the first term on the right has the form of a mean-field approximation. The second term can be expanded to 
give, up to second order in gradients of the order parameter 

PF[n]-pF{n^) J d^i Apf {^.-r {^1)) (23) 

-^yc2 (ri2;no) (pi (7*1; T (7*0) + pi) ^^^^^£^^^7*21 • Vr (7*1) ^7*1^7*2 

-I J C2 (ri2;no) (pi (7*i;r(7*i)) +pi) ^^^^'^^^^"^^^^ 217*21 • V VF (7*i) cj7*id7*2 

y C2 (ri2; no) (pi (7*1; r (7*1)) + Pi) ^'^gr2'|yj'^^ (^21 • Vr (7*1)) ' d7*id7*2 
OH use a specific parameterization 



pi (7*; r (7*)) = po ( 1 + To (7*) + 5^ exp (iKn ■ 7*) r„ (7*) j 

V n>0 ) 



(24) 



where Kn is the nth reciprocal lattice vector. Substituting into the expression for the free energy gives 

pF\ny\-pF{^y) ^ yd7*i A/?/(7*i;r(7*i)) (25) 



-^Po y C2 (ri2;po) (po + n^^o (7*i) -I- ni) 7*2i • VTo (7*i) dl^id~ri 
-^Po^] /'c2 (ri2;po) exp • 7*12) r„(7*i)7*2i • Vr„ (7*i) d7*id7*2 

-^Po y C2 (ri2;po) (po + PoFo (7*1) +pi)^2i7*2i • V VTo (7*i) (i7*i(i7*2 
-\p^^Y^ / C2 (r-i2;po) exp (^iiir„ • 7*12 j r„ (7*i) 7*2i7*2i • VVr„ (7*i) d7*id7*: 

n>0"' 
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The linear terms vanish by symmetry of the uniform fluid and, after changing variables and integrating by parts, the 
final result is 



f3F [m] - (3F (ni) ~ J dT^i A/3/ {T.-T (r* i)) 



+ ^Po I C2 (ri2;po) (~r*2i • To (7i))^ 



+ 



^Po [ C2 (?'i2;Po)exp (j,Kn ■ 7*12) (7*21 • Vr„ (7*1)) d'r'id 



r 12- 



(26) 



This can be written in the form 



PF [m] - I3F (ni) fdr Apf ; T (t)) + ^ ^ i^g' {OiT, (7^)) {dj^ )) 

L ab 



(27) 



with 



Kf-" = Sab^pl J exp (iKa ■ 7*) rirjC2 (r; no) d'r'. 



(28) 



The main criticisms of this theory are that (a) all correlation functions are eventually expressed in terms of an 
expansion about the uniform liquid and these expansions truncated at the two-body dcf; (b) the latter stages of the 
derivation depend on a particular parameterization of the density and (c) the mean field contribution, A/ (7*i; T (7*i)), 
is not a truly long-wavelength approximation as it clearly exhibits variations on the length scale of the lattice. To 
clarify the last point: in a long- wavelength theory, one would expect to see the term 



/ 



dti A/3/ (7*1; r (7*1)) 



replaced by something like 



j dR j dti A/3/ (7*1; r p)) ^ j P)) 



(29) 



(30) 



so that the variational equations for the order parameters do not contain terms that vary over microscopic length 
scales. 



III. EFFECTIVE GINZBURG-LANDAU MODEL 



Ideally, one would like to obtain a Helmholtz free energy in which microscopic length scales do not occur of the 
generic form 



fo (r (y)) + lifg^ (r (r)) f Ar, (r)) (^n i^) 



dr-i 



(31) 



in which dependence on the atomic structure has been "integrated out". Furthermore, one would expect the first 
term to take the typical mean-field form of the bulk free energy evaluated for the local order parameters, 



F/o (r (7^1)) = pF Hlr=r(7^.) = PFbuik (r (7^)) . 



(32) 



where (i^bulk (I^o) is the free energy of a bulk system with spatially uniform order parameters To. Given a free energy 
functional of this form, and truncating at second order in the gradients, the Euler-Lagrange equations, eqs.(13), would 
become 



92 



< (r (^)) ^r^r, (7^) + (r (T^)) ( ir^b (7) ) ( i^^. (7) 



dvidrj 



d_ 



d 



dr-j 



d 



dTa (T) 



/3Ffa,ifc(r(7*))+M<5ao = 0, (33) 
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where it is assumed that Tq is the average density. These variational equations are derived under the usual assumption 
that the values of the order parameters on the boundaries are held fixed. 

The idea presented here is to first break up the volume of the systems into cells containing a single atom. Then, 
withing the cells, the spatial dependence of the slowly varying parameters is expressed as a gradient expansion. Since 
these parameters are, a priori assumed to vary little over atomic length scales, it is assumed that the gradient expansion 
can be truncated without serious approximation. Finally, the resulting expressions for each cell are resummed thus 
separating the atomic and structural length scales. 



A. Initial gradient expansion 



To begin, it is convenient to use the representation 

pF[p]= j dr*/3/(T^;[p]) (34) 

where / i/r; [F]) can be expressed in terms of the one-body dcf and the density using the methods of the previous 
Section, 



/?/ (^; W) = - ci (^5 [\p]) p dx. 

Jo 



(35) 



An important assumption made throughout this Section is that the density is parameterized by a set of order pa- 
rameters which might include e.g. the average density but that the underlying lattice structure, including the lattice 
constant, is held fixed. For example, the Haymet-Oxtoby parameterization given in eq.(24) is of this form since the 
reciprocal lattice vectors Kn are held fixed; another example is that used in ref.([7]) 

P (^) = Vo (^) (^) ^ E (^) - ') (36) 

where the two order parameters are Fq (T*) = rjo (T*) and Fi (T*) = a (T*) but the lattice vectors -R„ are held 
constant. Furthermore, attention here will be restricted to the most important case of simple Bravis lattices, such as 
FCC and BCC, with a single atom per unit cell. Notice that this also includes as a special case the nonuniform liquid 
which results in both parameterizations from taking F„ ("r*) = for n > 0. These models are therefore sufficiently 
flexible to describe the solid-liquid interface, as well as e.g. the gas-liquid interface. 

To begin, partition space into the Wigner-Seitz cells where the center of the n th cell, Vn, will be at and all 
cells have volume V^s. The expression for the free energy can therefore be written as 



Nv n 

l3F\p] = Y, / drUfir-M) (37) 



Nv 

In each cell, the density is expanded about that of a system with uniform order parameter 

Nv 



0PiP] = E / dt pf{r;[pn]) (38) 

n=l -IVr. 

V / dr dn (5/?(ri;F(ri ,F„ 

(5p( r i;F( r i) ,F„)<5p( r 2;F( r 2) ,F„) 



n=l 



(5p(r 1) 



n=l 



5p{r i)5p{r 2) 



Pn 



+... 



where 



5p (y 1; F (T^i) , F„) = p (7^1; F (T^i)) - p (y 1; F„) (39) 
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and r„ = r is the value of order parameter at the center of the nth cell. At this point, all terms are retained 

and no truncation is imposed. Next, the dependence of the densities on the order parameters is expanded using 



p(ri;r(yi)) = p(yi;r„) + ^(r„(yi)-r„,a) 



dp {'r\;Tn) 



(40) 



ao 



+... 



as is the spatial dependence of the order parameters themselves 

r„ (^i) = r„ (tr.) + (t^i - Rn) ■ (vr,,)^ + (^^ - • ("^^^O^ + 



(41) 



The key assumption made is that terms of order V"r may be neglected for some n > hq. Then, the functional 
expansion, eq.(38), is only needed up to order no as well. Taking no = 2 and combining eqs. (37-41) gives 



Nv ^ 



(42) 



n=0 

Nv 



W(^;[r]) 



+ 



1^ /• ,^ /• W(^;[r]) 



dp{'ri;Vn) 



n=0 

TV,/ 



5p{r i) 



d^p{rv,Vn) 



+ ^ E / dl^ [ d'ridl^; 



W(^;[r]) 



5p{r {)6p{r 2) 



(Tx - Rn) (j2 - Rn) : ( Vr„,a) ( Vr„,b) 

ap(^i:r,,) ap(r2;r„) 



Because the terms involving functional derivatives are being evaluated at constant values of the order parameters, 
they exhibit the symmetry of a bulk solid. In particular, they are invariant with respect to translations by a lattice 
vector so a change of variables yields 



Nv . 

f3F[T] = V / dl' Pf(y;[pn]) 



(43) 



Nv ^ , 

+E/ 


dl^i 




1"' 


n=0 •' 




+^e/ 





W(^;[r]) 



(5p(ri) 

W(V;[r]) 



r 1 • vr„,a 



Op( T*i;r„) 



(5p(r 1) 

W(^;[r]) 



9p (T*i;r„) 



5p ( r 



r 1 r 1 • 



W(r;[r]) 



5p{r i)5p{r 2) 



dp (7*1; r„) 5yo(7*2;r„) 



The second term on the right vanishes by virtue of reflection symmetry in a Bravais lattice, p (— T^i; r„) = p {'r i; r„). 
The only explicit dependence on is via (T*; [p„]) and the integral of this quantity and its derivatives over Vq can 
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be extended to an integral over all space and using 

J Vn "J 



(44) 



L 



d r 



S/3f{r;[T]) 



d r 



5p{r \) 



5p{r i)5p{r 2) 



^ l_ dlSFe^ [p] 
N Sp{l^i) 

1 S^0F,, [p] 



N Sp{ri)5p{r2) 



-j^C2 {'ri,~r2; [pn]) 



gives 



, Nv 



(45) 



n=0 

Nv 



dp (7'i;r„) 

aV(^i;r„) 



ri=0 ' ' 

~2]V 2^ / 1^ '' 2 C2 ( r 1, r 2; [p„]) r 1 r 2 • (^Vr„_a J (^Vr„,| 



This represents the desired coarse-graining of the dependence of the order parameter on the microscopic structure. 
This is made clearer by writing the free energy as 



Nv 
n=0 



Ko (r ( i?„) ) + Ea "^i.a (r ( i?„) ) : ( V V r„) 



with 



i?2,ab(r(i?„)) = 



(46) 



(47) 



1 

"27V 
1 



d'ri ci ( 7*1; [p„]) 7*17*1 



dT^i ci (7*1; [p„])7*i7* 



r idr 2 C2 ( r 1, r 2; [Pn]) r 1 r 2 



dp (7*1 ;r„) 

g'pCF*i;r„) 

5p(7*i;r„) 5p(7*2;r„) 



B. Resummation 



Having obtained the form given in eq.(46), the question is whether the lattice sums might be written as inte- 
grals so as to obtain something like the postulated form given in eq.(31). This question has not been addressed 

in previous derivations. To do the resummation, we first write the lattice sums more explicitly. Let the primi- 
tive lattice vectors for the I?-dimensional lattice be {^i}^i so that the set of lattice vectors can be written as 

{ Rn} = \Y^i=i'm'i i( where the various limits, Mj and iVj, will depend on the geometry. 

L J n=l L J rni...mD=Mi...MD 

Then, for any function of position / (7*) one has that 

No Ni(m2...mri) 

E/P»)" £ - E f{mi^i + ... + mD^D), (48) 

n mo=MD mi=Mi(rn2...mD) 

where, as indicated, the limits of the inner sums must be allowed to depend on the indices of the outer sums so as 
to allow for the most general boundary conditions. For simplicity, the dependence of the limits on the indices will 



9 



be suppressed but it should always be considered to be present. This summation can be expressed as an integral by 
means of the Euler-Maclaurin summation formula giving the exact relation 



Ni 



/(mi'a 1 + ... + mr('aD) 



(49) 



mi=M: 



Ml 



fix'ai + ... + mD'aD) 



(r + l)! 



Br+l (x) 



■E 

fc=i 



(k + iy. 



^/ (x a 1 + ... + mn a d) 



QT+l 

Ml 



fix'ai + ... + mD'aD) 



dx 



where Bk is the k-th Bernoulli number, Bk (x) is the k-th Bernoulli periodic function and r is an arbitrarily chosen 
positive integer. The Bernoulli functions are polynomials for x G [0, 1] and are defined outside this interval by 
Bk (x) = Bk (xmodl) . Thus, they are bounded for all values of k and in fact it is also true that their integral over 
the unit interval vanishes. The idea here is that since only terms up to order no are retained, and since the functions 
Bk (x) are bounded, the last term on the right can be neglected provided r + 1 > no- The second term on the right 
is evaluated on the boundary. For r = 1 , it vanishes provided 



—f {x'ai + ... + ruD'aD) 



x=Ni 



d ^ ^ 
= g^f (x a 1 + ... + ruD a d) 



(50) 



x=Mi 



which could be true either because of the boundary conditions or because of symmetry. Even if it does not vanish, its 
contribution will be of order compared to the first term , where L = \ {Ni — M\) "a*i| so that in the thermodynamic 
limit, it should not contribute. Choosing r to be arbitrarily large and demanding that 



E 

fc=i 



{-if+^B, 



k+l 



{k + l)\ 



Qk 



fix'ai + ... +mD'aD) 



= 



Ml 



imposes a more complex boundary condition. To illustrate, in one dimension, r = 3 gives 

,3 



|i [/' (Niai) - f (Miai)] - ^ [/(^) (Niai) - /(^) (Miai) 



= 0. 



(51) 



(52) 



It could be argued that in the present case, since only terms up to second order in the gradients of the order parameters 

are kept, the third order terms in this boundary condition could in any case be neglected so that any choice of r > 2 
would be satisfactory. It seems best, however, to note that ,63 = so that taking r = 2 avoids any subtleties that 
might arise on the boundary. Similarly, if one were keeping gradients up to order no, and no is even, one would 
take r = no with the consequent boundary conditions involving derives up to order no — 1. No such simple choice is 
available if no is odd. 

So, for the present case of a second order theory taking r = 2 gives, when third derivatives are neglected, 



Ni 



/(mi'a 1 + ... + mr('a£)) ~ / / (x'ai + ... + mn 0*13) dx. 



mi=Mi 



Ml 



(53) 



together with the boundary condition, eq.(50). Furthermore, it happens that |i?3 (x)| < ~ 5 x 10~^ and that 

S3 (x) is oscillatory over the distance of one lattice spacing whereas the order parameters are meant to vary slowly 
over such small distances, thus further strengthening the argument for neglecting the higher-order term. When this 
is extended to include the second sum, a new complication arises since this gives 



^2 Ni JV2 pNi 

E E f {'mi'ai + ... + mD'aD) ^ / dy dx f [x'ai + y'a2 + ■■■ + niD'a d 

™ J Mo J Ml 



(54) 



m2=M2 mi=Mi 



El 
' 2 

1 

3! 



d_ 

dy J Ml 

N2 g3 

dyBa (y) — 

M2 dy^ 



dx f {x~ai + y'a 2... + niD'a d) 



N2 



M2 



Ni 



dx f {x'ai + y'a 2 + ••• + rno'a d) 



Ml 



10 



In the event that Mi and Ni depend on y, as they would e.g. for a spherical boundary, then the second term on the 
right becomes 



dx 



Ml 



d_ 

dy 



f {x'ai + y'a 2... + niD'a d) 



N2 



M2 



(55) 



+ 



dNi jy) 
dy 



f {Nilxi + y~a2— + niD'a d) 



dMi 
dy 



f (Mi'a 1 + y"a 2... + niD'a d) 



N2 



M2 



Now, the first contribution can be required to vanish giving a simple and expected multidimensional generalization 
of eq.(50) while the second term involves only the value of / on the boundary and geometric functions concerning 
the definition of the boundary. All of these quantities are fixed in the variational procedure leading from eq.(31) to 
the Euler-Lagrange equations, eq.(33), so they are of no importance. Thus, with the necessary generalization of the 
boundary condition, the second term on the right of eq.(54) can be neglected. Finally, the third term on the right in 
eq.(54) becomes 



I-N2 

JM2 



dyB3 (y) 



— jg-/ {Ni a I + y a 2 + - + niD a d) a 1 + y a 2 



dy^ 



niD'a d) 



(56) 



+3 



N2 



+ 



dy"^ dy 
'dNi 
dy dy"^ 



f {Nilti + 2/7? 2 + ••• + niD'a d) 



dyB3 (y) 

M2 

fN2 

+3 / dyBs (y) 

JM2 

dyBs {y) / dx —rf {x'ai + y'a 2 + ... + niD'a d) 

M, J Ml 



d'^Mi d 
dy"^ dy 



f (Mi'a I + y~a2 + ■■■ + niD~a d) 



dM d^ 

/(TViQ* 1 +y'a2 + ... + mD^D) ^^-^ (Mi'a 1 + y'a2 + ...+ niD'a d) 



Here, the first term is of no importance since it involves the value of / on the boundary and the fourth term can be 
neglected as it is of third order. Note that in the remaining terms, ^/ {Ni'ui + y~a2 + ... + niD'a d) is a derivative 
evaluated for variations on the boundary, as are all of the other derivatives that occur in the second and third terms. 
Thus, if the boundary condition is that the order parameters assume constant values on the boundary, then these 
terms are negligible. This type of boundary condition is certainly reasonable for spherically symmetric problems or 
on the radial boundary for those with cylindrical symmetry. On the other hand, for planar interfaces one could take 
the volume to be defined by Ni = —Mi = N for some constant N in which case = 0, etc. and again the second 
and third contributions are negligible. Again, these terms should be of no consequence in the thermodynamic limit 
as they are of order compared to the first term in eq.(54). Thus, with these considerations, it has been proven that 
up to third order gradients, 



N2 Ni 
m2=M2 mi=Mi 



mD~a d) 



N2 



dy 



M2 



Ni 



dx f {x'ai + y~a2 + ■■■ + mD~a d) ■ 



Ml 



Clearly, this result can be extended to any number of dimensions to get the general result 



niD 



'aD)= / ••• / dxi...dxD f {xi'ai + xd'S d) 

J Mn J Ml 



(57) 



(58) 



/ Md J Ml 

which is expected to hold as the limits go to infinity or, for finite systems, with the corresponding boundary conditions 



5 _> 
Vz, — / (xi a 1 + ... + XD a d) 



XD~a d) 



(59) 



and / is cither a constant on the boundaries or the boundaries are defined by fixed values of Mi and Ni. Since the 
primitive lattice vectors are linearly independent, the integral can be written as a volume integral by a change of 
variables giving 



No Ni 

- E f{mi^^ + -+^D-^D) = ^ [f{l')+0{V^f)]dr + x, 



(60) 



mD=MD mi=Ni 



where x represents the neglected terms which involve the value of / on the boundaries and the various geometric 
factors. It is important in terms of calculating the absolute value of the free energy, but plays no role in formulating 
the Euler-Lagrange equations for the order parameters. 
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C. Continuum Limit 



Making use of eq.(60) and NVyjs = V, one has that, up to second order in the gradients and an overall additive 
constant, that 



PF [T] 



dR PF 



W j 



(Pr idPr 2d R 012(1^1 ■,^2; p(r(^iijj j 



dR ar„ (e) 



r 1 r 2 : 

dR dR 



dVai R] 



dVj R] 



An integration by parts gives 

/3F[r] = 1 1 dR(iF[p{T{R)) 



W j 



d'r^id'r^2dR 



C12 [ r I, r 2; 



p(r(i?jjjj (^1^2 - T^iT^i) : 



ar„ R) 



/ d^^dR ^ . c, (r,; [p (r (i?))]) 



'5r„ (r^ dp(rv,r(R)y 



dR 



or, rearranging the second term and using Gauss's theorem on the third term, 

1 , / /_>x X dTa (r) dVb (r 



f3F 



[r] = / 



dR 



1 



dR,, dR, 



1 

'2V 



1 • 



'dTa (r) dp(rr,T(^R 



S{V) 



dR 



with 



^ijO^) = ^/ d'r'-^d't'i ri2iri2jC2 (T*i,T*2;pr) 



dTa (r) 

dp{ri;T) dp{r2;T) 



■ dS 



dTa 



dTi 



(61) 



(62) 



(63) 



(64) 



(65) 



and where the surface term involves an integral of R over the boundary of the integration volume V , S{V). Except for 
the surface term, this is the desired result. In the thermodynamic limit, the surface term gives negligible contribution, 
as long as the order parameters and their derivatives are finite on the boundary. 



IV. CONCLUSIONS 



In this paper, it has been shown that, up to second order in the gradients of the order parameters and neglecting 
various constants which depend on the value of the density on the boundaries of the region of interest, the Helmholtz 
free energy of an arbitrary simple crystalline system with fixed lattice structure can be written in the form of eq.(31) 
thus generalizing the results of Lowen, et al[5],[6]. This expression makes a cleaner separations of length scales than 
does the older form, eq.(27). Furthermore, at no point in this derivation were higher order correlations arbitrarily 
neglected: the only approximation made was to truncate the gradient expansion at second order. This directly 
addresses some criticisms of the use of these types of models, namely that they are too crude in the treatment of 
correlations[10]. 
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The applicability of cq.(64) depends on the boundary conditions. It requires no qualification in the thermodynamic 
limit, except that it is assumed that the order parameters and their derivatives are finite on the boundary and the 
usual derivation of the Euler-Lagrange equations require that the values of the order parameters be constant on the 
boundary. For finite systems, the resummation of the local free energies requires that 



d ^ ^ 
Vz, — r (cci a 1 + ... + XDaD) 



Xi=Ni 



d ^ 
— r (xi a 1 + ... + XDa d) 

OXi 



(66) 



Xi=Mi 



and that F is either a constant on the boundaries or that §^ 

OXi 



, etc. Problems involving spherical or cylindrical 
symmetry will generally automatically ensure the constancy of F on the boundary and planar interfaces can be handled 
by assuming a volume defined by Ni = —Mi = N for some constant N. For finite systems, the surface term in eq.(64) 
must also be considered. 

One noteworthy point about the Ginzburg-Landau form, eq.(64), is that all of the microscopic quantities are evalu- 
ated for spatially uniform order parameters or, in other words, for bulk systems. Thus, subject to the approximation 
inherent in the gradient expansion, this technique provides a framework for using good models of bulk systems to 
study interfacial systems. 

One important problem not addressed here is that of elastic relaxation. In general, one would like to allow the 

lattice parameter to be included in the list of order parameters so that the lattice can expand or contract as one 
passes from, say, a bulk solid phase into a bulk liquid phase. In fact, this model has been used for a similar problem 
in which the lattice is allowed to vary from FCC to BCC[11]. The present derivation is not sufficiently general to 
allow for this as it would mean allowing the volumes y„ to depend on the order parameters. This should be possible 
but the results can be expected to be considerably more complex than those presented here. Such a generalization 
will be the subject of a future publication. 
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